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Abstract 

We study an axially symmetric solution of a vortex in the Abelian- 
Higgs model at critical coupling in detail. Here we propose a new idea 
for a perturbative expansion of a solution, where the winding number of 
a vortex is naturally extended to be a real number and the solution is 
expanded with respect to it around its origin. We test this idea on three 
typical constants contained in the solution and confirm that this expansion 
works well with the help of the Pade approximation. For instance, we 
analytically reproduce the value of the scalar charge of the vortex with an 
error of O(10 -6 ). This expansion is also powerful even for large winding 
numbers. 
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1 Introduction 


A significant feature of many gauge theories is the existence of topological soli- 
tons which may appear when the gauge and/or global symmetries are spon¬ 
taneously broken. Monopoles, vortices and domain walls are by now familiar, 
and have found important applications in vast areas of modern physics, such as 
cosmology, condensed matter physics and particle physics. From another view 
point, the topological solitons can be seen as nontrivial solutions of nonlinear 
differential equations. A direct way to study topological solitons is solving such 
nonlinear equations exactly. For instance, a beautiful systematic method to 
construct exact solutions for instantons has been well-established and is widely 
known as the ADHM construction (T. This is, however, a special case and for 
many other types of solitons numerical calculations are needed to study solu¬ 
tions. 

The present work concerns the so-called Abrikosov-Nielsen-Olesen (ANO) 
vortex as the simplest topological soliton with finite energy in the (1 + 2)- 
dimensional theory. This vortex appears as a topological defect {3] in Ginzburg- 
Landau theory [2] and may be viewed as a static solution to the equations 
describing the 1+2 dimensional Abelian Higgs model [Jj. In this theory all the 
vortex features depend on one dimensionless parameter A = m s /mj [ ]: the ratio 
of the Higgs boson mass m s to the vector boson mass m v . The intervortex 
force is, roughly speaking, a superposition of an attractive force caused by the 
Higgs boson and a repulsive force caused by the vector boson as seen in a scalar 
potential m 

U{R) - 7^ {~qlK 0 {m s R) + q 2 m K 0 {m v R)) (1.1) 

for a well-separated pair of vortices with a large distance R. Here, q s and q m 
stand for a vortex scalar charge and a magnetic dipole moment, respectively. 
Therefore the force with the longest correlation length is dominant and the vor¬ 
tices attract (repel) each other for A < 1 (A > 1) [6]. The critical coupling A = 1 
is a rather special case where net intervortex forces are exactly canceled thanks 
to the coincidence of the two coefficients, q s = q m = 27tCi. From a mathe¬ 
matical viewpoint, the Euler-Lagrange equations reduce to the first order dif¬ 
ferential equation called the Bogomornyi-Prasad-Sommerfield (BPS) equations 
for vortices saturating Bogomol’nyi bound, whose total energy is quantized as 
Ek = \k\TTV 2 with the winding number k £ Z. In this critical case, the constant 
Ci appears, for instance, in a potential for a pair of moving vortices m 

U\=i(R,u) ~ ttv 2 x ClK 0 (m v R)\u\ 2 + C(|u| 4 ) (1.2) 

with a relative velocity u 1 since only the magnetic field accepts a Lorentz boost 
and the two forces are not canceled out. Unlike the remarkable cases of in¬ 
stantons and monopoles, no analytic solutions for this BPS equation in flat 
spacetime have been found even at this critical coupling. Thus only a few quan¬ 
tities are exactly calculable and a detailed study of the vortices, for instance, 
the calculation of a value of Ci requires numerical analysis. 

In this paper, to complement the numerical analysis, we propose a simple and 
straightforward, but new idea for analyzing vortices at critical coupling, where 

1 A/\/2 is known as the Ginzburg-Landau parameter. 
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fields are expanded perturbatively with respect to the winding number k € Z 
around its origin k = 0. To justify this perturbative expansion, (let us call 
it “small winding-number expansion ”), the BPS equations must be extended 
so that they allow a real winding number k G R. Since the BPS equations 
with an infinitesimal winding number \k\ -C 1 can be exactly solved, we can 
systematically perform perturbation calculations without tuning any parameters 
and this perturbative expansion is supposed to work well as a practical tool. 
Here, we calculate values of three typical quantities with A = 1 including C\ as 
the most simple examples to check this idea. 

The constant Cj has often been calculated in the literature. De Vega & Scha- 
posnik [5] gave a semi-analytical study for axially-symmetric solutions with an 
arbitrary winding number k € Z>o, and constructed power-series expansions 
around a center of a vortex and asymptotic expressions for the opposite side. 
These two can be determined by only one constant for the power-series ex¬ 
pansion and Ck (Zk in their notation) for the asymptotic expression. Comparing 
these parameters in a middle region, they obtained the values: C\ = 1.7079... 
and D\ = 0.72791.... These values now seems to be widely accepted in lit¬ 
erature, for instance, C\ = 1.7079 appears in Refs. mini nans] and also in a 
standard textbook of Vilenkin & Shellard (8]. However, we encounter a different 
value for C\\ C\ = 10.58/27T ~ 10.57/27T ~ 1.682 ~ 1.684 which was obtained 
by Speight [TO] about twenty years later than de Vega & Scliaposnik [5j. Fur¬ 
thermore, Tong m gave the supergravity prediction C\ = 8 1 / 4 ~ 1.68179... 
which seems to agree well with Speight’s C \. These values also seem to be 
accepted in literature, for instance, Ref.[T2] and another standard textbook by 
Manton & Sutcliffe jj]. There exists a 1.5 % discrepancy between old and new 
results. 

In Sec.2.6, we shall conclude that the correct value is the old one C\ = 1.7079 
by using two different kinds of numerical calculations with higher accuracy. In 
Sec.4.2.3, we reproduce this value by using the small winding-number expansion 
to verify its power. 

This paper is organized as follows. In Sec.2, we review the BPS vortex in 
the Abelian-Higgs theory, and define an extended vortex function which allows 
the winding number of non-integer, as a solution of the BPS equations. There 
a non-trivial integral formula including the vortex function is derived and three 
typical constants C v , D v and S v for a vortex solution are introduced and their 
analytical and numerical properties are discussed. In Sec.3 we perform a small 
winding-number expansion of the vortex function and the three constants using 
Feynman-like diagrams. Results obtained there are modified in Sec.4, using 
the Pade approximation to overcome problems with finite convergent radii of 
the expansions. Summary and discussion are given in Sec.5, and some useful 
inequalities and details of the calculations are summarized in the Appendices. 


2 


2 Review of ANO vortex at critical coupling 

2.1 Set up for ANO vortex 

The Abrikosov-Nielsen-Olesen (ANO) vortex is an elementary topological soli- 
ton in the 2+1 dimensional Abelian-Higgs model 

c = ~h F ^ F>lv + t 2 - 1 ) 

where </> is a complex scalar field, metric is rj = diag.(+l, —1, —1) and covariant 
derivative is + i . A scalar potential V{(f>) is of the wine-bottle type 

VW = ^fm 2 -v 2 )\ ( 2 - 2 ) 

which has a vacuum \<p\ = v where the U{ 1) gauge symmetry is spontaneously 
broken. The Higgs mechanism makes the scalar and the gauge fields massive. 
Their masses are given by, m s = A ev, m v = ev respectively. The spontaneously 
broken U{ 1) symmetry gives rise to a soliton which is topologically stable object 
supported by tti(U (1)), of which element is called a winding number. To require 
vanishing of the kinetic term \T>i(j)\ 2 = 0 at the spatial infinity connects this 
winding number with the first Chern class 

7Ti([/(l))=Z 9 k = — J d 2 xFi 2 . (2.3) 

This topological defects are called the Abrikosov-Nielsen-Olesen vortex. 

In this paper, we take the critical coupling constant, A = 1, as the simplest 
model, where the two masses are identical, m v = m s = m. Then we can perform 
the Bogomol’nyi completion of an energy density Ji for static configurations as 

^L=1 = ^{A 2 ±y(t; 2 -M 2 )} + 1 -\{V 1 ±W 2 )^\ 2 

+y+ 12 ± '-Pdi {4 >v 3 4>) , (2.4) 

and a total mass (tension in higher dimension) of vortices, T, has a lower bound 

T= J d 2 xH\ x=1 > +y J d 2 xF 12 = ±nv 2 k. (2.5) 

The inequality is saturated by BPS states which satisfy the BPS equations 

+ F 12 = y {v 2 - \ct>\ 2 ) , (Pi ± iP 2 )0 = 0. (2.6) 

Without loss of generality we will consider the BPS equations with the upper 
sign. In order to find general solutions of the BPS equations, it is useful to solve 
the second equation in Eq. (12+1) at first and, it can be solved with the complex 
coordinate z = xi+ix 2 and introducing a smooth real function ip Teg = ij; ieg (z, z) 

k 

Az = -dzip ieg , <j) = ve ~P(z), P{z)=J\{z-zi), (2.7) 

7=1 
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where an arbitrary holomorphic function P(z) can be set to be a monic poly¬ 
nomial without loss of generality. Here zeros {zi = x) + ixj £ C} of the Higgs 
field (p are topological defects and identified as positions of vortices. One of the 
important features of BPS vortices is that they feel no interactions since the 
attractive and repulsive force are exactly canceled. So we can put BPS vortices 
anywhere as many as we like. Note that the smooth field ip Teg must behave as 
Vkeg ~ log | P(z)\ 2 at the spatial infinity to obtain a finite energy, it is convenient 
and more familiar to rewrite ip ieg in terms of a singular field ip, 


V’ = V’reg-log|P(^)| 2 = -log^—, (2.8) 

tr 

so that ip vanishes at the spatial infinity. With this singular field, then, the first 
equation in Eq. (EH) can be rewritten to be, so called, Taubes’ equation 

— + m 2 (l — e - ^) = J, (2.9) 

with source terms J 


k 

J = J(x) = dTr^Or- xi). 
1=1 


( 2 . 10 ) 


Here we used that the magnetic field can be rewritten as, 

- F 12 = 2 d z d s ip Teg = i ( dip) + J) (2.11) 

which coincides with Eg. (12.3D and k is the total winding number. Existence 
and uniqueness of a solution for Taubes’ equation with a given arbitrary J have 
been established by [14] . With this solution, therefore, we obtain a complete 
solution for <p and A,. In terms of a solution of ip and the source J, the energy 
density 'Kbps for BPS vortices can be rewritten to 

Heps = ^ (J + dfa[tp]) , a[tp] = ip + e - ^ - 1 > 0, (2.12) 

which gives the lower bound in Ea. (l2.5D . There is, however, no known exact 
solution for this equation, even in the simplest case with k = 1. 


2.2 Extension of Taubes’ equation and particle description 

In a case that kj vortices coincide at x = Xi for each /, the source terms are 
replaced with 

J = 47t^^ kjS 2 (x — Xi), k = ^^kj. (2.13) 

i i 

where Ap indicates the winding number at x = xj. A request that the winding 
number ki is positive integer is to give the single-valued Higgs field cp and Profiles 
of ip and the magnetic field in Eq. (12.IIP and the energy density in Eq. (12.121) can 
be calculated without constructing cp. If we omit constructing <p, therefore, we 
can formally extend Taubes’ equation with the generalized source terms 

J = 47r^ x/), vi £ {v \v > — 1, v £ R}. (2-14) 

i 


4 






Here the winding number kj is renamed iq to stress that vj can be non-integer 
and the lower bound of the winding numbers will be discussed in Sec.2.4. A 
‘total mass’ of this extended object is formally calculated as 

?bps = f d 2 x'H.BPS = nv 2 x i/, u = ^ zq, (2-15) 

•’ i 

which takes a negative value for v < 0. Integrating the both sides of Taubes’ 
equation Eg. (12.91) we find the following identity corresponding to Eg. (12.31) 

v = ^ J d 2 x (d 2 if + J) = r ^~ J d 2 x{ 1 - e - ^), (2.16) 

which is no longer an element of 7 Ti(/ 7(1)). In the rest of this paper, we will study 
this extended Taubes’ equation with a generalized source term Eo. (l2.14|) and 
its solution numerically and analytically. This extension allows us to consider a 
Taylor expansion of the solution with respect to the winding numbers as discussed 
in Sec(3] although we are not specially interested in the solution with the winding 
numbers of non-integer. 

Uniqueness of the solution for this extended Taubes’ equation can be easily 
shown as Appendix IA.1I For instance, we know the trivial solution 

0 = 0 for J = 0. (2.17) 

To show existence of the solution for the extended Taubes’ equation is difficult 
and out of scope of this paper, and we just assume the existence of the solution 
here. Therefore, the solution of if is a function with respect to a coordinate x, 
positions of vortices {x/} and their winding number {iq}, 0 = 0(x, {x/, zq}). 
Furthermore we assume that the solution is differentiable with respect to {zq}. 
Under this assumption, we can derive, for each I, 

(— d 2 + m 2 e~^)^~ = 47t<5 2 (x — x/) (2-18) 

ozq 

from Taubes equation Eo. (12.91) with the source Eq. (I 2.14|) . According to Ap¬ 
pendix IA.1I the above equation show that the solution 0 is strictly increasing 
with respect to each zq, dip/dvi > 0. In the limit of the vanishing source J = 0, 
furthermore we find 

lim ^ — -S 2 (x-xi) = 2A' 0 (m|x - x/|), (2.19) 

J->o ovi —of + m. 

where the modified Bessel function of the second kind Kq(x) emerges as a two- 
dimensional Green’s function. That is, in this limit a vortex solution is exactly 
solved and treated as a linear combination of free massive particles and for small 
\uj\ <C 1 at least, 0 is approximated well everywhere as 

0 « 2 ;//A'o(?n|x — x/|). (2.20) 

/ 

This is the starting point of the small winding-number expansion which will be 
discussed in SecEJ 
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In this particle description, it will be convenient to rewrite Taubes’ equation 


as 


00 ( _ 1 \n 

- d 2 3 ip + m 2 ip = J + m 2 a[ip], a [ip] = ip 2 ^ (2.21) 

(n + 2)! 


n—0 


with cr[/0] as dimensionless self-interaction terms. Then, by applying the Green’s 
function method to Taubes’ equation, we obtain^ an integral equation for ip with 
Green’s function G(x) = K 0 (m\x\), 

^(x) = 2'^2u I G(x-x I ) + m 2 J ^-G(x - y)a[ip(y)]. (2.22) 


Since a[ip] > 0 and K 0 (x) > 0 are always hold, we find that the solution of 
Taubes’ equation must satisfy a fundamental inequality 

ip(x) > uiKo(m\x — xi\). (2.23) 

/ 


2.3 Scaling argument and a physical size of a vortex 

Let us consider the following Lagrangian in a two-dimensional Euclidean space- 
time 

£bps = - \{diipf - m 2 {ip + e~^ - 1) + Ji/>, (2.24) 

which induces Taubes’ equation as an equation of motion of ip, and an actiorjf] 
is 


K = - 


d x£bps 


solution 


T -^ghost? 


(2.26) 


where iLghost is introduced to cancel UV divergences of the kinetic term and the 
source term and we set ILghost as, for instance, 

d 2 x Q(0i<p) 2 + ^m$<p 2 - Jpj , 

ip(x)=2'^2viK 0 (m 0 \x-xi\). (2.27) 

I 

After this regularization we can apply the scaling argument to this action. 
For simplicity, let us consider an axially symmetric case with the source J = 

2 Here we used the fact that t/j vanishes at the spatial infinity. 

3 Substituting the solution, K becomes a function with respect to complex coordinates 

zj = describing positions of vortices. With a limit of mo —> 0, this quantity gives a 

Kahler potential describing the vortex moduli space m as 

7rv 2 vi\zi\ 2 +v 2 lim K. (2.25) 

mp—> 0 

At the limit mo —»• 0, 1/mo gives a IR cut-off and -Rghost can be eliminated by Kahler 
transformations. Actually one can confirm that the above Kahler potential gives Samols’ 
metric 1171. 
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47 tv8 2 {x). Since K is a dimensionless quantity, the dimensional argument tells 
us 


n 2 dK 
0 = mr ——77 + 


dK 

dm 2 ' ‘ "°dm% 


m n 


(2.28) 


By using equations of motion for ip and ip, derivatives of K with respect to 
masses can be calculated by 


TO 


dK 
dm 2 
2 dK 
^dri 2 


J d 2 xm 2 [ip + e ^ — f ) = m 2 J cC 2 xip — 47 r v. 


2 ^-^-ghost 

l ° dm 2 


= - d 2 X^ 2 = - 27 TV 2 


(2.29) 


where we used Ea. (l2.16l) . Therefore, we find the following formula [T 8 ] 

/ 2t r 

d 2 xip = —77 x v(y + 2). (2.30) 

m 2 


As we seen the above this exact formula does not come from topological ar¬ 
gument, but from the scaling argument. To check numerical calculations we 
use this formula in this paper. Thanks to this non-trivial identity combining 
Ea. (l2.16)l . the following integral is calculated as 


/ 


~Hbps 

7TV 2 


1 

87T 


/ 


d 2 x\x\ 2 dfa[ip] = 


2tt 


{ip- 


o -V> 


- 1 ) = 


9 ’ 


(2.31) 


and a size of the vortex with the positive winding number v > 0 can be naturally 
defined with the energy density Kbps given in Ea. (12.12l) and calculated as, 


-Kbps = 



/ (Px\x\ 2, Hbps 
f d 2 x’H.BPS 



for v > 0 , 


(2.32) 


which turns out to be a key point in Sec HI It is natural for the scaling argument 
to determine a typical size of a soliton. 


2.4 Axially symmetric solution 

Let us consider a single vortex sitting the origin with the winding number p, that 
is, we consider a solution with the source term J = 47tp<5 2 (x). Its configuration 
is axially symmetric and described by a function ip = ip(mr 1 v ) with respect to 
a radial coordinate r = |x| and the winding number v. The partial differential 
equation m, therefore, reduces to an ordinary differential equation 


1 d 


dip 


rfrVTr^^ 1 -^ 


with the following two boundary conditions 


ii», M = - 2 ., 

r—>-0 dr 


lim 'ip = 0 . 

r—too 


(2.33) 


(2.34) 
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Even for the non-integer number v, a set of the differential equation and the 
boundary conditions defines an unique solution under assumption of its exis¬ 
tence. Especially for small \v\ -C 1, if is approximated in the full range of 
r € R>o as 

if ss Ei[ijj\ = 2uKo(mr), lim — = lim ^ = ifi = 2K 0 (mr). (2.35) 

v — ^0 V OV 

See Fig JT] for some examples of profile functions of N[if] which denotes if calcu¬ 
lated numerically. Here we assume that the solution if is smooth with respect 



Figure 1: Magnetic flux in the left panel and differences Ai/j = N[ip] — E\ [?/.'] in the right 
panel 


to v at v = 0. This assumption requires the solution to be extended for the 
negative winding number v. Since dif/dv > 0 as discussed in Sec l2.21 a lower 
bound of v is shown by taking a derivative of the both sides of Eq. (12.301) as 

0 < [ < e x W = ±L W+l) , (2.36) 

J ov m* 


that is, there exist no solution of Taubes’ equation with v < —1. We just assume 
the existence of the solution with v > — 1 in this paper. 

Note that we can show the following inequalities although we have no exact 
solution. Applying the discussion in Annendix lA.il to Taubes’ equation with 
the source J in Eq. ()2.14p . we find the solution if must be positive for u > 0 
and be negative for v < 0, and Eg. (12.331) tells us that is strictly increasing 
(decreasing) with respect to r for v > 0 (v < 0), and therefore the boundary 
conditions Eo. (12.34l) give lower and upper bounds as, 


if > 0, —2v < r^y- < 0, 

dr 

dif 


for v > 0, 


if < 0, — 2v > r— >0, for — 1 < v < 0. 

dr 

According to Appendix IA.11 the following inequality 

( ^ +ra v*)g_ m .^(£)*>o. 


(2.37) 


(2.38) 


implies that if is a downward-convex function, 
























Combining Eq. fl2.35l) with this fact, we find that ij)/v is strictly increasing with 
respect to v and furthermore we obtain 

dib ib 

—— > — > 2 Ko(mr) >0 for v > 0, 

ov v 

0 < ^ < — < 2K 0 (mr) for — 1 < v < 0. (2.40) 

ov v 

With this axially-symmetric solution ip(x) = V’( r ) with r = |ir|, the integral 
equation Ea. (12.22l) reduces to 


pOO 

ip(r) = 2vKo(mr) + m 2 / ds sGp(r, s)cr[+s)], 

J o 

where the reduced Green’s function Gp(r, s) takes the following form 


(2.41) 


Gf(t,s) = J —Kq ytnsjr 2 + s 2 — 2 rs cos 6j 

= 0(r — s)Ko(mr)Io{ms) + 0(s — r)Ko(ms)Io(mr) (2.42) 

with the step function 0(a;) and the modified Bessel function of the first kind 
I 0 (x). 


2.5 Observable parameters, C v , D u , S u 

2.5.1 Dv and Internal size R m 

To define the solution if} of Taubes’ equation even with the positive non-integer 
winding number v, we have to consider a behavior of the solution around the 
core of the vortex seriously. Note that in the massless limit m —> 0, Taubes’ 
equation has a general solutior@ with a positive real arbitrary constant R[ n as, 


lim x!) = - log Y, Y = ( — ) (2.43) 

m->0 V Irin J 


and with the finite mass m > 0, therefore, ij) can be expanded by m and we find 
an expansion of ij) around the origin r = 0 in an unfamiliar form, 


V’ 


-log Y +Y, F n(Y){mr) 2n 

n= 1 


—2^ log (mr) + 2Djy + 


4 

-2D 


4 (1+1/ 


k ( mr ) 2 

?z(mr) 2 ( 1+I ' 


for v > 0 
for - 1 < v < 0 


(2.44) 


where we treated mr and Y as if they were independent of each other, and a 
function F n (Y) is independent of m and turns out to be a polynomial of order 
n with respect to Y determined sequentially by solving Taubes’ equation as, 


F " {Y) 4 i 1 (1 + +)’ 64 ((2 + v) 2 (1 + +) (2 ' 45) 

4 Here we omit the boundary condition for the spatial infinity. 
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which must vanish in the limit v —> 0 for a finite radius r due to Eg. (12.171) . The 
dimensionless constant D v appeared in the expansion is related to R m a^H 


D v = vlog(mR in ). 


(2.47) 


Therefore the expansion of i /> can be defined by a pair of parameters {v,Ri n }. 
The uniqueness of the solution with a given v means, however, that to satisfy 
the boundary condition at the spatial infinity, the constant R m must take a 
certain value corresponding to each value of v, that is, a function i?j n = Ri n (p), 
otherwise a function defined by the expansion always glows up at a large r. In 
Appendix IA.21 this feature is analytically discussed and at the present we find a 
pair of lower and upper bounds of R m as 


2y/^ + T 

m 


> R\n 



for v> 0 . 


(2.48) 


According to Ea. (l2.40|) R ln and D v /v turn out to be strictly increasing functions 
with respect to v and take values at v = 0 


lim —- = lini(Ar 0 (^) + logr-) = log2 — 7 « 0.115932, 

v —^0 V r —>-0 


lim R ln = 

v —>o 


2e " 7 1.12292 


(2.49) 


m m 

with Euler’s gamma 7 . In Fig[2] we plot a profile of D v /v. Note that there is 




Figure 2: Profile of D v for the full range of v. Numerical Data ;Y s hi [tA] are plotted by 
dots. Dashed lines in the left panels describe E n [D v ] given in Sec [71 Dashed lines in the right 
panel give the bounds given in Eg. 112. 48 ll . t ’2 [ D ;/ ]. P(, \D U \ plotted by a solid line are defined 
in SecQ] 

an another way to calculate D u using the integral form Ea. (12 .41(1 as, 

fib \ m 2 f°° 

D u = lim (-hi / log(mr) I =i/(log2 — 7)4 -/ dssA' 0 (ms)cr['!/'(s)](2.50) 

r-s-o y 2 ) 2 J 0 

5 A relation between D v for v = k S Z>o and defined by de Vega & Schaposnik [5j 

is 

jeM-ZDk), (2.46) 

For instance, we numerically obtain D\ = 2 exp(—2 X 0.505360825 ...) = 0.72791247.. . which 
coincides with their value Df = 0.72791. 
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These different two definitions of D „ will be used to double-check numerical 
calculations of D v . 

Since the axially symmetric vortex solution we consider has the only one 
mass parameter to, we expect that the dimensionfull parameter R ln controlling 
a profile of the solution should be the same order of the vortex size -Rbps given in 
Eq. (12.3211 . Thanks to Eq. (12.4811 , roughly speaking, we find actually i?BPS ~ R'm 
for large v. We call an internal size. On the other hand D v is directly related 
to a value of the action K with J = 47 tvS 2 ( x ) in the previous subsection. In the 
same way of Eq. (12. 2911 . we can calculate a derivative of K with respective to v, 

= — 47 t x lim (ip — ip) = 8i tv log (] — 877 ( D v — p(log 2 — 7 )) (2.51) 
dv r->- 0 \ m o / 

and by setting the mass of the ghost too to be too = 2 e _ 7 m, we obtain the 
following simple relations, 


D„ = - 


1 dK 
87 t dv ’ 


K = -8 tt 


dyD v . 


(2.52) 


2.5.2 Scalar charge C„ 

Let us take to large conversely, that is, consider a infrared region r ^ R ln ss 
2 yjv/m. There, an asymptotic behavior of 7 /’ can be treated as a fluctuation of 
a free massive scalar field around the vacuum. Due to the axial symmetry, such 
a fluctuation is written with a certain constant C v € R>o as 


ip S 3 2C u Ko(mr). 


(2.53) 


There is the similarity between this asymptotic form and the form of Eo. (12.3511 
and the uniqueness of the solution of Taubes’ equation indicates that the two 
constants C„ and v are in one-to-one correspondence. Actually, to satisfy the 
boundary condition at the origin r = 0, the constant C v must be a function 
with respect to v and according to Eg. (12.171) . Eg (12.191) and Eo (12.391) we find 


lim C), = 0, 


.. dC v 

lim —— 

v-¥ 0 dv 


= 1 , 


d 2 C v 

dv 2 


> 0 . 


(2.54) 


These property tell us that C v /v is strictly increasing with respect to v and a 
lower bound of C v is given as C„ > v. A profile of this function is shown in 
FigEl According to the integral equation Eo. (12.41l) . C v can be calculated by 


c v = lim ——j —r 
r->-oo 2Ko(rnr) 


to 2 

= v -\—— / dssIo(ms)<j[ip(s)]. (2.55) 

2 Jo 


Bringing this identity back, we can remove the explicit p-dependence from the 
integral equation Eo. (12.41D as 


pOO 

ip{x) = 2C v Ko(mr) - / ds sG^r, s)a[ip(s)] : 

Jo 

with an ‘advanced’ Green’s functiord 

G a d(n s) = 0(s - r) {K 0 (mr)I 0 (ms) - I 0 {rnr)Ko{ms)} > 0. 


(2.56) 


(2.57) 


6 Positivity of this quantity is easily shown since Ko(r) ( Iq(t )) is strictly decreasing (in¬ 
creasing) with respect to r. 
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Figure 3: Profile of C u for small v in the left panel and for large v in the right panel. 
Numerical Data N 8 bt[Cj/] are plotted by dots. Dashed lines in the both panels describe 
approximants of the order n, E n [C u \, in terms of the winding-number expansion discussed 
in Sec|3] Ps[C u ] plotted by a solid line and Pi[Cu] plotted by a dot-dash-line are defined in 
Sec [4] 


Using this integral equation Ea. (12.56l) . the asymptotic behavior in Eg. (12.531) is 
modified as 

/•OO 

'ip = 2C v K 0 (mr)-2Cl / ds sG ad (r, s)K 0 (ms) 2 + 0(e~ 3mr ). (2.58) 

Jo 

Thanks to these two different forms of the integral equations for ip Ea. (12.41[) 
and Ea. (12.56D . we find lower and upper bounds as 

2uKo(mr) < ip < 2C v K d {mr). (2.59) 

A one of purposes of this paper is to confirm the true value of C\. 

2.5.3 Total scalar potential S v 

Finally let us consider the following definite integral 




(2.60) 


which is dimensionless and proportional to a total potential energy of the 
Abelian-Higgs model at critical coupling, 


S v = 


A dE„ 


E 1 d\ 


A=1 Ei 


d 2 xV{(p) 


A=l,sol 


(2.61) 


This quantity with v > 0 satisfies 


n „ ni / ax. 

0 < < Ty 2^r (1_e ) = I/> 

and according to Eg. (12.171) and Ea. (l2.19D we find 

lim S„ = lim —A = 0, lim " = 2. 
v -¥o v-¥0 dv i/->o dv z 


(2.62) 


(2.63) 


7 This quantity also appeared as a fundamental constant, c = 2Si ~ 0.830707,in Eq.(5.2) 
of a paper ESI- 
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Thanks to Ea. (l2.40|) we find that S v is also an increasing function with respect 
to v and according to the profile of S v shown in Fig0] an ‘energy’ per an unit 
winding number S^/v is also an increasing function with respect to v, and this 
property gives 


S„ 1+V2 > S Vl + S V2 . (2.64) 

This inequality is consistent with the well known property of type H (type I) 
vortices, that is, intervortex forces are repulsive (attractive) for the coupling 

A > 1(A < 10. 


S(v) SM 




Figure 4: Profile of Su for the full range of v. Numerical Data N r s ht[<S'i/] are plotted by dots. 
Dashed lines in the left panels describe E n [Sv] given in Sec|3| Ps[S u ] plotted by a solid line 
and Pi [S u ] plotted by a dashed-line are defined in Sec(4] 


2.6 Numerical Data 

We numerically calculate values of G V ,D V1 S V in most of the range of v as 
v = 1/20,1/10, • • •, 500,1000 using mainly the shooting method. These data are 
listed in Table 03 We will denote these data as Wht [C/], iVsht [£0] and -/V s ht[S/] 
for C V ,D U ,S U respectively. In SecU we use these data as references to show 
how the winding-number expansion introduced in Sec|3] works well. The other 
purpose of this subsection is to settle the problem on the numerical value of C\. 
We need, therefore, numerical calculations with high accuracy. To show accu¬ 
racy of our numerical data to readers, let us enter into details of the numerical 
calculations we performed. 

Note that there exist two kinds of strategies in the shooting method and 
we observe a big difference in usability between them. We calculate numerical 
solutions of '0 in a region {r|e < r < L} where we set m = 1 and take e = 
10 -2 ” +1 and L = 2 v / i/+plogl0 with p 1 n = 8 ~ 9 referring to the flux size 
-Rflux given in Ea. (12.32D . The first strategy is to take r = e as the initial 
point of the calculation and fine-tune the parameter D v so that a profile of 
0 satisfies the boundary condition at r = L and read C v from a profile of 0 

8 It is natural to expect the following inequalities on values of total energies Ej. for axially- 
symmetric vortex-solutions, 

Ek 1 +k 2 < + E k2 for A < 1, (2.65) 

which induces the inequality ( 12.6211 . To the best of our knowledge, there is no known mathe¬ 
matical proof for these inequalities although they are quite reasonable. 
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V 

c v 

D u 

Rin 

Su 

1/20 

0.05152300 

0.007221252 

1.155375 

0.002320344* 

1/10 

0.1061386 

0.01714170 

1.186986 

0.008668711* 

1/5 

0.2249350 

0.04429221 

1.247899 

0.03070642* 

1/2 

0.6633334 

0.1736933 

1.415364 

0.1444002* 

1 

1.707864 

0.5053608 

1.657584 

0.4153533 

2 

5.336582 

1.443305 

2.057831 

1.085081 

3 

11.86421 

2.615596 

2.391367 

1.832041 

4 

22.61080 

3.948209 

2.683313 

2.619544 

5 

39.31961 

5.402536 

2.946174 

3.432922 

10 

317.5504 

13.88300 

4.008030 

7.704638 

20 

5424.053 

34.27687 

5.550253 

16.68079 

50 

1284274. 

107.9305 

8.659094 

44.65765 

100 

5.455139x10 s 

250.0538 

12.18905 

92.38242 

200 

2.607156xl0 12 

568.9475 

17.19704 

189.1678 

500 

4.568733xl0 19 

1650.717 

27.15154 

482.7929 

1000 

6.065189xl0 27 

3647.519 

38.37932 

975.6104 


Table 1: Numerical Data of C v , D v {) and S„. All data are sufficiently stable values and 
we double-checked them except for data added stars. 


at r = L. Since the initial conditions are given by a pair {is, D v }, an incorrect 
pair always makes a profile function blow up at large r. The second one is 
to take r = L as the initial point and fine-tune the parameter C v so that 
v = —(rip'/ 2) at r = e and read D v at r = e. In this strategy the profile 
function is controlled by the only one initial parameter C v which is related to 
v in one-to-one correspondence thanks to dC v /dv > 1. With the sufficiently 
large L, therefore, a profile function with an arbitrary C v always gives a certain 
solution corresponding to a certain v , without the profile blowing up, and thus 
this strategy gives a function v = f(C v ). Thanks to this property, it is easy 
to create a computer program for tuning C v automatically with a given is and 
arbitrary precision. We take the second strategy in this paper although the first 
strategy was taker@ in Speight’s paper US¬ 
As we explained above, numerical data N s \ lt [C v ], Wht [D v \ for C V ,D V are 
directly obtained. To double-check those data, we also use the integral for¬ 
mulas Eq. (12.5511 and Eq. (12.5011 for C v , D v respectively, to obtain different data 
KulCvlKulD.}- We regard |Ad ht [X]/Wht[A'] - 1| with X = C„,D„, as errors 
of these data and plot them in the right panel of Fig|5] For instance, we obtain 
as double-checked numbers, 

AshtlCi] = 1.707864175 

7V s ht[£>i] = 0.505360825378 ( 2 . 66 ) 

for is = 1 and the numerical data listed in Tabled] have been double-checked in 
this sense. Therefore we conclude that the numerical result C\ = 1.7079 given 
by de Vega and Schaposnik is correct. Thanks to the non-trivial identity in 

9 He stated there as “Hence, all numerical solutions blow up at large r, and even though aj 
and 62 were tuned to six decimal places, the Runge-Kutta algorithm could not shoot beyond 
r = 10.” 
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Eq. (I2.30[i , we can estimate accuracy of the profile functions itself by calculating 
the following quantity 


5 = 


1 


v(v + 2) 


drrN[ip] + 2N[C V ] j dr?'iv 0 (r)| — 1 


(2.67) 


and we plotted this in the right panel of Fig[5] Note that we observe that the 
precision of jV s ht[Cy] generally get worse than those of S,D U as shown in Fig|5] 
The precision of calculations in Speight’s paper seems to be less than six digits 
and we guess that his result C\ ~ 1.683 has an error of 0(1O -2 ) ~ 0(1O -3 ) 
which is consistent with the other numerical results including ours. 


i Error of C 
Error of D 
. Error of S 




nr 9 



• Error of C 
■ Error of D 

♦ 6 




Figure 5: Estimated numerical errors: The left panel plots errors of numerical data calculated 
by the relaxation method from those calculated by the shooting method as, |A r r i x [X]/jV s ht[A] — 
1| with X = Ch, D v . S„. The right panel plots errors of numerical data in terms of the shooting 
method itself as, |N' ht [C^]/iV s ht[C',,] — 1|, |A r ' ht [T ) i.]/A r s ht[^ ) i/] - 1| and S. 


We obtain also a stable numerical value of Si with long digits 

Na,t[Si] = 0 . 4153533072562 , ( 2 . 68 ) 

by the shooting method. To perform double check of the values of S„, we also 
use the relaxation method as the other numerical calculation. In the relaxation 
method, we introduce a relaxation time r and extend ip(x) to be dependent on 
r, -0 = ip{x,r), and modify the equations of motion by adding a friction term 
dip/dr with an appropriate signature. With an appropriate initial function 
of ip, ip(r,T = 0) = 2vKq(t) for instance, this friction term defines the time 
evolution of ip and decreases an ‘energy’ of this system defined in Eg. (12.261) . 
In principle, therefore, the true solution could be obtained with an infinite r 
as tp(r) = lim T _ >00 ip(r, r). As larger r, we will get better accuracy in many 
cases. In reality, beyond a certain finite r, we observe stability of values of the 
observables with small noises, since those accuracy can not be better than the 
calculation accuracy. For instance we stopped the time evolutions at r « 4 x 
10 4 . The relaxation method is convenient and powerful to solve (simultaneous) 
nonlinear (partial) differential equations numerically. We need no fine-tuning 
of any parameters there. In the simple system we are considering, however, 
the shooting method is more powerful to get precision. Generally speaking, 
numerical data iV r ix[A’] for X = C V ,D V , S v calculated by the relaxation method 
get worse precision as shown as Fig J5] We find |A r r ix[<S’i]/A r s ht[«5'i] —1| ~ 5xl0 -11 
which is guessed to be mainly an error of iV r i x [S'i]. We also get iVrixfCi] = 
1.707864188 ... and Wix[-Di] = 0.5053608253753... again. 
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3 Small Winding-Number Expansion 

In the paper J5}, de Vega and Schaposnik calculated C\ and Di by a semi- 
analytical study. Their strategy was essentially as follows. Let us divide the 
integrals in Ea. (12.501) and Ea. (12.55l) as 

nb pOO 

+ , with b ss i? B PS ■ (3.1) 

Jo Jb 

The former integral is calculated by inserting the expansion Ea. (l2.44|) which 
depends on D u and the latter is calculated by the expansion Ea. (12.581) which 
depends on C„. Then we obtain simultaneous equations for C„ and D v , and 
thus, approximate the values of C Vl D v as their solution. 

In this section we will give a different expansion of the solution ip using 
Eq. (12.411) and calculate them more straightforwardly and more systematically. 

3.1 z^-expansion of the vortex function ip 

In the normal case, we can not define an expansion of ip with respect to the 
winding number as a topological quantum number. In the previous section, 
we relax the winding number v from an integer to a real number and assume 
smoothness at v = 0, and thus, we can consider a Taylor expansion of the 
solution for ip with respect to the winding number as, with Ea. (l2.17D 

OO 

ip =^2v n ip n . (3.2) 

n= 1 

Since the approximate solution E\[ip] = uipi in Eq. (12.351) satisfies the boundary 
conditions Ea. (12.341) and has the same asymptotic form as Ea. (12.53l) for an 
arbitrary z/, we expect that the following finite series of order n 

n 

E n [iP}=^v m ip m (3.3) 

m—1 

gives a good approximation and becomes better as the larger order n. Here, a 
higher-order coefficient ip n for n > 2 can be sequentially calculated by expanding 
the integral equation in Ea. (12.22l) . or Ea. (12.411) for the axially symmetric case, 
with the first approximant E\ [ip] , as 

/ j2 oo 

-^G(x - y)cr n (y), a[ip] = ^ v n a n (3.4) 

n—2 

where expansion coefficients o n = cr n (x) in the interaction terms cr^] are 

^2 = ^1. cr 3 = -^1 +^1^2, •••• (3.5) 

Let us call this Taylor expansion a small winding-number expansion , or simply, 
a ^-expansion. Note that in this expansion the winding number v is fixed and 
higher order corrections have no logarithmic singularity as 

lim = 0 for n > 2. (3.6) 

i —s-o dr 
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The absence of the solution for v < — 1 shown in Eg. d2.36p might indicate that 
a radius of convergence for the ^-expansion of ij) is less than 1. In Sec HI we will 
discuss that this fact is not a big problem. 

We can perform calculations of the ^-expansion of i/> with the familiar technic 
using Feynman diagrams. The ^-expansion of i/j(x) is given concretely as 


ip(x) = 2u x x 


0 +i(2^) 2 x 



o 


0 


+(2 vf x 





+eV), (3.7) 


using conventions for Feynman diagrams, 


G(x) = K 0 (m\x\) = x 


'I ^- G ^-y) G (y) 2 = 



Here diagrams of the order n have n external legs coming from the point-like 
vortex at the origin x = 0. 


3.2 E n [C v \ 

Let us approximate C v analytically by using the i/-expansion, 

OO 

C„ = J2 c n" n , a = I- (3-9) 

n= 1 

In principle, its coefficients c n can be obtained by taking the ^-expansions of the 
both sides of Eg. (12. 5511 and inserting i/j n obtained in Ea. d3.7l) into the right hand 
side. Comparing Ea. d2.dll! and Eo. d2.55p . however, we find that the coefficient 
c n can be calculated by only replacing the propagator with Io(m\x\) as 


fpn{x) = 




(3.10) 


.(3.8) 
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where the triangle symbol stands for 


I 0 {m\x\)= x -. (3.11) 

For instance, coefficients 02,03 are calculated as 


C2 = 


C3 = 


See Appendix iBl for details. Finally we obtain 

2 

C v = u +-^=u 2 +—u 3 + 0.0126799A 1 
3 a/3 108 

—0.0013557(41)i/ 5 + 0.000781(22)i/ 6 + 0(v 7 ), (3.13) 

which gives a finite series E n [C„] as an approximant of order n 

n 

E n [C v \=Y J CkV k - (3.14) 

fe =i 

As shown in Fig| 6 ] we observe that as the order n is larger, an error of E n \C\], 
that is, |.Era[Ci]/-Afsht[Ci] — 1| is smaller. The sixth order approximant for v = 1, 
Eq[C]\, gives a quite nice value near to the numerical value fV s ht[Ci] in Ea. (12.6611 
as 



Eq[C{\ = 1.70809..., 


£ 6 [Ci] - iVUtCr] 
IVshtlCr] 


1.0 x 10“ 4 . 


(3.15) 


Unfortunately the accuracy of this value is worse than that of the value Ci ss 
1.7079 given by de Vega and Schaposnik. According to FigH>]a radius of conver¬ 
gence of the infinite series, u c , is obviously finite and smaller than ten, v c < 10 
and we can not judge whether v c is larger than one or not. In Sec]4j we will 
overcome these problems. 


3.3 E n [D v \ 

Next, let us consider the ^-expansion of D„, 

OO 

D v = ^2 dnV n , d\ = log 2 — y. (3.16) 

71=1 
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Figure 6: Errors of the n-th order approximants E n \C \] and P n [C\) in the left 
panel, and E n [C\o] and P n [C\o\ in the right panel. Pn[C u \ will be defined in 
Sec0] 


According to Ea. (12.5011 . the expansion coefficient d n for n > 2, is calculated by 
reducing diagrams in Eq. lETTb as, 


4>n(x) = 




(3.17) 


We find therefore, by performing integrals numerically, 



= 0.115932z/+ 0.585977i/ 2 - 0.333905+ 3 + 0.244999// 

—0.196695// + 0.165065(79)// + (D(v 7 ) (3.18) 


and the //-expansion of R ln is also obtained as 

r r> \2 (‘2Dv\ 

( mR in ) = exp I —— 1 

= 1.26095+ 1.47777//+ 0.0238675// 2 

-0.030728// + 0.0300632// 4 - 0.02652(10)// + 0(// 6 ). (3.19) 

Note that this quantity is known to have the lower bound 4e -1 // and the second 
coefficient is near to this bound as 1.47777 > 4e _1 = 1.47152. Finite series 

n 

E n [D„}=Y J d k v k (3.20) 

k=1 

were expected to be good approximations, but we find their slow convergence 
as seen in FigH 
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3.4 The ^-expansion of the formula Eg. (12.301) 

To check consistency of the {/-expansion of the formula Eq. (12.301) , we need some 
unfamiliar formulas. There is a non-trivial identity as, 


J d-xi/j n — 

and using Eq. JBU) we find 


cPxm 2 
—d 2 + m 2 


a n = d 2 xa ni 


(3.21) 


crxi/j-ii/jn = / d 2 xa n 


m 


- d 2 


-ipi = — m 2 / d 2 x<j, 


dipi 

dm 2 


. (3.22) 


Using the above formula, we also find with o\ = 4tt 6 2 (x)/m 2 


/ ^2 = \ [ d 2 x= -2 tt|^ 


= ^ (3.23) 

r=0 m z 


and since is a dimensionless quantity we can confirm 

1 


d xi/j 3 = / d~x -~4>i + V’l^ = d x [ --i/h - 0^1 m 


6 


1 . 3 1 2 _2 ^4 


6 


1 9 


6 9m 2 


m 2 / d 2 xipf ) = 0. 


dm 2 / 
(3.24) 


To check Eo. (12.301) for more higher order, similarly we must need the dimensional 
argument again. Checking Ea. (12.30l) is, therefore, tautological in this sense. 


3.5 E n [S v ] 

To calculate the {/-expansion of S v , at first we rewrite the definition of S v by 
inserting the identity in Ea. d2.31l) 


S' = ^/^(l-e-) 2 = ^/^(, 2 -,3 


+ •••) 


v 2 + m 2 f ^ Q(l- e -^) 2 + l_ e -V’_ V , 


d 2 x ( V’ 3 V’ 4 7 , K 1 


{/■“ + 77U / -I — ——b —-U 0 H-U 0 - 

' 2 tt V 3 4 60 ^ 24 ^ 2520 


1 P 7 + 0{'i/; 8 ) . (3.25) 


Here we canceled a 1 /j 2 term to avoid complicated and redundant calculations 
such as those in Sec l3.4l and thus, substituting Eq. 113771) we easily find the fol¬ 
lowing expansion^ 

OO 

Sv = y^s k v k 

k=l 

0 0 

= v 2 — -(2{/) 3 x + (2{/) 4 x 

o 



10 Here a diagram of order n has n + 1 external legs. 
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+0(u 6 ) (3.26) 

and then, we obtain by reusing the calculations of integrals in Ea. (l3.18l) 

S v = v 2 - 1.562605;/ 3 + 2.73802i/ 4 

—5.05307iz 5 + 9.59699 p 6 - 18.5461(5)i/ 7 + 0(u 8 ). (3.27) 


A finite series of order n for S v is defined as 

n 

E n [S v ]=v^s k v k , S! = l. (3.28) 

fc=l 

Unfortunately we find, however, that these finite series do not work as approx¬ 
imations even at v = 1 as shown in Fig0] and it is inevitable to use some 
technique for obtaining good approximations. 


4 Pade approximations and Large v behaviors 

4.1 The bag model for large v 

The result of the vortex size Tbps i n Eg. (12. 321) implies that the total magnetic 
flux of a vortex is proportional to an area occupied by the flux for v > 0, 


/ 


d 2 xFi2 


lib 2 

2nu = — x ttR bps 


( 4 . 1 ) 


where m 2 /2 = e 2 v 2 /2 is the maximum of the magnetic field allowed by the 
BPS equations Eg. (12.61) for v > 0. This fact evokes the liquid droplet model of 
nuclear structure, and gives an intuitive explanation in our axially symmetric 
case for the Bradlow bound m, which means just that the area 7 t1?bps mus t 
be less than the total area if we considered a closed two-dimensional base space. 

In a paper |23] , the size 7 ?bps w &s obtained by a physically intuitive way 
using the bag model proposed in EHH2 for the large winding number u. In 
the bag model, a vortex configuration consists of an inside Coulomb phase, 
the outside vacuum in the Higgs phase, and a thin domain-wall at r = R 
interpolating their phases. In the Coulomb phase, the magnetic field takes a 
non-vanishing constant determined by the total magnetic flux in Eq. m with 
v = k, and vanishes in the vacuum. By omitting a thickness of the domain-wall, 
profiles of the Higgs field and the magnetic fields are approximated by 


I0I 2 


0 for r < R 
v 2 for r > R 


1*121 


for r < R 
0 for r > R ’ 


of which the total energy is calculated as 
2’kv 2 


2bag — 


2 4 


e 2 R 2 


nR 2 > 7 tv 2 x v = Tbps- 


(4.2) 


(4.3) 
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Figure 7: Configurations of the magnetic flux = 1 — e ^ for u = 9,36,81,144 of 

which radiuses are estimated to be m.R = 6,12,18, 24 respectively. 


This energy is minimized just at R 2 = Av/e 2 v 2 = -Rj|pg- Actually, we numeri¬ 
cally observe profiles of the magnetic field for large v in FigJT] A profile of the 
domain-wall is almost invariant with various values of v. For large v, therefore, 
a contribution to the total energy Tb ag form the domain-wall can be negligible. 

Since a vortex configuration for large v drastically changes around the domain- 
wall at r ~ R 1/to, we expect that the approximation for r -C R ln in 
Eg. (12.431) is applicable for r = R — e < R with e = 0(l/m) as 




—2v log (m(R — e)) + 2D„ 



(4.4) 


and similarly the asymptotic behavior in Eq. (12.531) is applicable for r = R + e 


0(1) « MR + e) » C„K 0 (m(R + e)) « (4.5) 

V 2rnR 

Inserting R = i?sps = 2y/v/m, these estimations give large-i/ behaviors of C v 
and D„ as 

C v « 0(1) x uie 2 ^, A,»^(logi/ + 0(1)). (4.6) 

We also estimate S v as 

S„ = I -^—-(1 — e~^) 2 ~ —— x 7T (mR) 2 + y- x 0(1) x 2irmR. + O(R 0 ) 

2 J 27t 47t 47t 

» V - PVv + O(v 0 ). (4.7) 

Not that the term proportional to \fv comes from contribution of surface of the 
vortex and the coefficient /3 must be positive due to Eg. (12.621) . The above esti¬ 
mations for large v will become important clues to modify the approximations 
using the ^-expansion. 


4.2 (Global) Pade approximations 

Let us assume that we know only a finite series of order n, 

n 

E n [F{v)\=Y J hv k , (4.8) 

k=0 
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as a part of a certain infinite series F(v) and it behaves as almost an alternating 
series like F{v) = |/o| — \fi\v + \f 2 \v 2 — ..., and it seems to have a small radius 
of convergence v ss u c . To get a good approximation for v > v c with such a 
series, it is powerful to use the Pade approximation which replace the series by 
some rational functions, with n = m + l, 


F n [F{y)] = P {m , l) [F(u)\+0(u m+l+1 ). (4.9) 

where a Pade approximant of F(v) is given by 

m 

ao + ^2 anl,n 

P { m,i)[F(u)} = --, (4.10) 

1 + J2 b rF n 

n =1 


where coefficients of these rational functions are determined so that they satisfies 


d k F(y) 


dv k 


„=0 " 


i /=0 


for A; = 0, ra + /. 


(4.11) 


Here the two sets {a n } and {&„} are determined uniquely from the finite set 

{/o; /21 ■ ■ ■ I fn+m}- 

There is arbitrariness in a choice of {rn, l) for the order n. The approximant 
P(m,i)[F{v)\ behaves for large v as 


P(m,l)[FM] « fV 


™ ,,m—l 


(4.12) 


Note that if we fix p = m — l to remove that arbitrariness, then n is restricted 
so that n — p = 21. In the case of p = 1 for example, we arrange the Pade 
approximants for all of the order n as 


P( 1 , 0 ) [W], 


yJP i2 ,0)[ F m P(2,l)[P(v)}, ( 4 - 13 ) 


4.2.1 P„[S„] 

The series expansion for S u seems to be almost alternative series, and according 
to configurations for the finite series E n [S u \ shown in the left panel of Fig0] 
we guess that the radius of convergence is around \u\ ss 0.5 which implies, for 
instance, that the function S„ has a singularity at v ~ —0.5. The Pade ap¬ 
proximation can avoid such a singularity and enlarge the radius of convergence. 
Let us take the following rational functions P n [S u \ with respect to v, as Pade 


approximants 

of the order n for 

P2[S V \ 

= P(2,l)[Sv] = 

P 3 [S„] 

= ^(3,1) ['S'i/] = 

Pa[S „] 

= -P(3,2) [^i/] = 

P 5 [S„] 

= P( 4,2)[S„] = 

P e [S„] 

= -P(4,3) [^i/] = 


v 


1 + 1.56261/ 
v 2 +0.189609j/ 3 
1 + 1.752211/ ’ 
v 2 + 1.05188i/ 3 
1 +2.61449i/+ 1.34739i/ 2 ’ 
v 2 + 1.34536i/ 3 + 0.0556454+ 1 
l + 2.90796i/+ 1.86162i/ 2 ’ 

v 2 + 1.94979j/ 3 + 0.69144J/ 4 
1 + 3.5124^ + 3.44191z/ 2 + 0.814411i/ 3 ' 


(4.14) 
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Here we have fixed arbitrariness on choice of the Pade approximants P( m . n ) 
so that all coefficients of the above are positively As a result poles and zeros 
of these functions turn out to sit only on the negative real axis of v as shown in 
Fig El and the rational functions P n [S v ] have poles around v ss —0.5 in common. 
Actually these functions give good approximations in a wider range of v as 


• pole of P[S] 
■ zero of P[S] 


• pole of P[D] 
■ zero of P[D] 


Figure 8 : Poles and Zeros for P n [5„]/z/ 2 in the left panel, and for P n [D u \/u in right panel. 
We observe common poles v « —0.5 for P n [S u ] and v « — 1 for P n [D] and a common zero 
v « —0.2 for P n [D\. A pair of a pole and an adjacent zero do not change a large-z^ behavior 
remarkably. 


shown in Fig (9] Note that these rational functions behave as 


Errors 



Figure 9: Profiles of P n [Su] in the left panel and their errors from numerical data of S u 
calculated by the shooting method, \Pn[Su]/N s ^[Si/] — 1| in the right panel. 


Pin[S v \= 0{v), P 2n +i[S v \ = 0(y 2 ), for large v, (4.15) 

and P 2 n[Sv] give comparatively good approximations even for large v. This 
property can be understood if we take account of the behavior of S v for large 
v shown in Eq. (14.71) . Extra zeros v « —5.27 of Pi [SV] and v « —23.4 of P 5 [S„} 
shown in FigJ5]can be regarded as disturbances for large- 1 / behaviors. 

Let us consider the large- 1 / behavior more seriously. The large-i' behavior 
in Eq. (1X71) does not always mean that the function S v has a branch cut. For 
an example, a function y / Ftanh(y / F) has no branch cut anywhere although it 

11 This fact might be just by our good luck. We have no proof for existence and uniqueness 
of such a choice in the all order n. At least, we have to avoid zeros and poles on the positive 
real axis of v since we know 0 < S v < although the arbitrariness remains under this 
restriction. 
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behaves yd/ for large v £ R>o- Here, we just assume existence of a branch cut. 
For instance, a function 


has a branch point at v = —1/2 and desirable behaviors as 

SroiJ - 2 + 0(- 3 ) for v « 1/2 

1 " l + X ) for i/» 1/2 


(4.16) 


(4.17) 


and consequently it works as a quite good approximation for the full range of v 
as shown in Fig|3] The Pade approximation taking account of informations for 
large v is called the global Pade approximation [24] . Note that an expansion of 
the following quantity is also alternative series due to the singularity, 

^1-— j = 1-2 I / + 4.12521 i/ 2 - 8.60125i/ 3 

+18.0239i/ 4 - 37.857i/ 5 + 79.5748 j/ 6 + 0(u 7 ), (4.18) 


Let us apply the Pade approximation to the above series or its squared quantity. 
According to Ea. (14.7D . the above quantity behaves as 0(i/ _1 ) for large v and 
this property fixes the arbitrariness of Pade approximants completely. Addition 
to Pi[S v ] in the above, then, we obtain the following functions as the global 
Pade approximants of S v , 


h[s v ] 


h 3 [^] 

P 4 [S„] 


h[Su} 


h[Su] 


i 


l + 4z/ + 3.74958z/ 2 ’ 


v — 


v — 


V 


V 


1 + 0.80192^ 

1 + 2.80192i/+ 1.47863i/ 2 ’ 

1 + 0.697034z/ 

1 + 4.69703 v + 6.53772i/ 2 + 2.31356^ ’ 


v — 


v — 


V 


V 


l + 1.11774i/+ 0.064997i/2 
1 + 3.11774i/ + 2.17527t/2 + 0.0904502 j/ 3 ’ 

1 + 1.81492i/ + 0.525555z/ 2 
1 + 5.81492^ + 11.5348J/ 2 + 8.60739i/ 3 + 1.63522+* ’ 

(4.19) 


which behave for large v as 


Pn[S v \=V-p n yfi + 0 



(4.20) 


with coefficients for n = 1, 2, • • 

{/3„} = {0.707107,0.718628,0.736437,0.740872,0.847699,0.75294, • • •}. (4.21) 

At this stage we do not know whether f3 n converges to a true value of /3. As we 
see in FiglTUl the global Pade approximation works well and Pe[S„] has a quite 
small errors less than 10~ 3 in the full range of v. Even for small i/, the global 
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Figure 10: Errors |Pn[S' ! ,]/fV s } 1 t[>S , id — 1| of Global Pade approximations P„ [.S',/ for S u . 
Distortion of a profile with n = 6 , at v = 1/20 is consistent to errors of A/^ lt [SV] itself shown 
in Fig [5] 


Pade approximants P n [S v \ give the best result as shown in FigfTIland the best 
approximant Pq[S v ] gives 


P 6 [5i] = 0.4153585..., 


Pe[5i]/lVsht[5i] - 1 


1.3 x 1CT 5 . (4.22) 


These are the satisfactory values enough as results with the small winding- 


Errors Errors 



Figure 11: Errors of S i in the left panel, and .S'-, in the right panel, 
number expansion. 


4.2.2 P n [D v ] 


The z'-expansion of D v given in Ea. (13.18l) also seems to be almost an alternating 
series and have a finite radius of convergence as shown in Fig[2] Hence let us 

12 We wish, although, to modify a slow convergence of the large-z/ behavior if possible. Note 
that a natural and probable expansion of S v around the infinity v = oc is 


Si, = v — /3y/v + ^2 

71 = 0 



(4.23) 


although our global Pade approximants P n [Su] set Q 2 n = 0 . If an actual expansion has non¬ 
vanishing Q! 2 n, convergence of P n [Su] is interfered by this feature. An irregular behavior of 
P 5 [(Sy] shown in Fig l 101 might be caused by this obstruction. This technical difficulty might 
be fatal unfortunately. 
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consider Pade approximations of E n [D v \. We can fix arbitrariness of the Pad’e 
approximation by requiring that all coefficients are positive as, 


Pz[D v ] 

Pa[D v ] 

Ps[D„] 

P e [D u ] 


0.115932!/+ 0.652038!/ 2 
1 + 0.569826!/ ’ 

0.1159321/+ 0.67104!/ 2 + 0.0960493!/ 3 
1 + 0.7337391/ ’ 

0.115932!/ + 0.706900i/ 2 + 0.297736i/ 3 
1 + 1.04306!/ + 0.176257!/ 2 
0.1159321/+ 0.728018!/ 2 + 0.419974i/ 3 + 0.0174966i/ 4 
1 + 1.22522!/+ 0.309917i/ 2 


which have a pole !/ ss — 1 in common as seen in FigJS] As shown in Fig |121 
P n [D v ] give comparatively good approximations. To get better approximations, 


D(v) 





Figure 12: Profiles and errors of P n [D v \ 


let us apply the Pade approximation not to D v it self, but to exp(2 nD v /v) with 
btain 

"0.853276 + v s 


n = 1, 2,4, then we 

Pi[D v ] 

= \ log 

h[D v \ 

= \ log 

Pi[D„\ 

= \ log 

h[o v \ 

= \ log 

h[p v \ 

= \ log 


0.676695 

V0.708551 + 1.66078!/+ !/ 2 ' 

0.667557 

0.654559+ 1.609821/+ !/ 2 
0.519101 + 0.668311!/ 

^0.511978 + 2.400061/ + 4.25790!/ 2 + 3.38280i/ 3 + i/ 4 

0.670835 

2.17939 + 5.34930!/ + 4.17049!/ 2 + i/ 3 


1.72838 + 2.21671!/ + 0.676831i/ 2 
These functions have the same behavior for large v as Eq. m, 


Pn[D„ 


■log 


a n < e/4 « 0.679570. 


(4.25) 


(4.26) 


Hence, as shown in Fig[l3]and FigflTl these give quite good approximations and 
errors of P$[D V \ are less than 10 -3 in the full range of v. The best approximant 
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Figure 13: Errors in global Pade approximations for D„. 


we obtained gives 

PelDi] = 0.5053639 ..., |P 6 [£>i]/AW-Di] - 1| » 6.1 x 10~ 6 , 

2 exp(-2P 6 [Di]) = 0.727908..., (4.27) 

which reproduces the numerical result presented by de Vega and Schaposnik. 
This value with the similar accuracy was also obtained analytically in Ref. j25|.@ 








Figure 14: Errors of D\ in the left panel, and D 5 in the right panel. 


4.2.3 P n [C v ) 

The ^-expansion of C v in Eg. (13.131) gives a quite good approximation for C v , at 
least for v = 1 and we do not know whether the radius of convergence is larger 
than one or not. In this stage, therefore, it is not useful to apply the (ordinary) 
Pade approximation to E n \C v ). Once we take account of the large-zz behavior 
of C„ given in Eg. (14.61) . however, we notice that there exists a singularity at the 
infinity and we have to remove this at the first stage. 

Let us consider the following function 

C v = sinh( 2 \/z/) (4.28) 

which has an infinite number of zeros on a negative real axis of v and regular 
everywhere except for an essential singularity at the infinity. The nearest next 
zero to the origin is v = — 7 r 2 /4 k, —2.47. It is, therefore, natural to assume 
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that a quantity F v = {C V /C U ) A has an infinite number of poles (and zeros) on 
the negative real axis of v. Actually we find that an expansion of F v gives an 
almost alternative series as, 


F v = 



1 - 0.248268^ + 0.020833^ 2 + 0.034017+ 3 
—0.0342630^ 4 + 0.0226871i' 5 + 0(v e ). (4.29) 


According to Eg. (14.61) . F v must behave for large v as 


F v 


const. 

v 


+ 0(v~ 2 ), 


(4.30) 


which means that we removed the singularity at the infinity in success. Next, 
let us apply the Pade approximation to the series in Eq . (14.29l) or its squared 
quantity, satisfying the property Ea. (14.30l) . We obtain! 13 ! 


Pi[C v \ 

= a 

P 3 [C„] 

= c. 

Pa[C v ] 

= a 

P 5 {C„\ 

= c. 

P e {C„] 

= c, 


P 2 \C V \=C„ 


1 

1 + 0.248268^’ 


1 

1 + 0.496535^ + 0.143244i/ 2 ’ 

1 + 0.712165^ 

1 + 0.960432^ + 0.217611iz 2 ’ 

1 + 0.600743^ 

1 + 1.09728^ + 0.441534^ 2 + 0.0481954^ 3 ’ 


1 + 0.709639^ + 0.0702914^ 2 
1 + 0.957906^ + 0.287275^ 2 + 0.017348+*' 


(4.32) 


where we added P\ [C u \ to the above although it does not satisfy Eq. (14.301) . We 
observe the large-iz behaviors of them except for P\[C V ] as 


Pn[C n ] 






(4.33) 


with coefficients 


{w 2 , w 3 , • • •} = {0.354169,0.318735,0.336252,0.342689,0.354693, • • •}. (4.34) 

In Figjl5j we observe that these functions give nice approximants in the full 
range of v and modify E n \C v \ as shown in FigJOj Resultantly, even for Ci, we 
succeeded in reproducing the numerical result C± = 1.7079 given by de Vega 
and Schaposnik as 


P 6 [Ci] = 1.7078629..., 


PejCi] 

iVehttCr] 


= 7.2 x 10~ 7 . 


(4.35) 


13 There exists still arbitrariness on a choice of a function C u . We can choose, for example, 

C v = z^cosh(2\/T). (4-31) 

However, a Pade approximant of the 6-th order with this choice turn out to brake up due to 
emergence of zeros or poles on the positive real axis of v. 
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Figure 15: Errors of global Pade approximants P n [C v \ for C v . 


5 Summary and Discussion 

We considered the small winding-number expansion (the ^-expansion) of the 
solution of the Taubes equation by extending the winding number, which is a 
topological quantum number, to be a real number larger than —1. We confirmed 
that the ^-expansion is useful to give good approximations of axially-symmetric 
vortex solutions in most of the range allowed for the winding number. Finally we 
found that for the scalar charge C\ the best approximate value in terms of the 
I'-expansion with the help of the Pade approximation is /^[Ci] = 1-7078629 ..., 
which coincides with a value A r s ht [Ci] = 1.707864175, obtained numerically by 
the shooting method. We judged that the result given by de Vega and Schapos- 
nik is correct, and Tong’s conjecture giving C\ = 8^ ~ 1.68 from superstring 
theory perspective is incorrect as a vortex solution in the Abelian-Higgs model. 
Their numerical similarity might suggest a certain universality. 

The Abelian-Higgs model of critical coupling is just the simplest toy model 
to test and establish usefulness of the ^-expansion. The idea of the j/-expansion 
is rather simple and more straightforward than the strategy taken by de Vega 

6 Schaposnik. As for BPS states of vortices in further complicated systems 
like non-Abelian gauge-Higgs models or of separated (parallel) multi-vortices, 
therefore, it is expected that the ^-expansion can be straightforwardly applied 
to their analytical approximations. Since it is difficult to apply the shooting 
method to such complicated systems, we guess that the role of the ^-expansion 
will become more important there. The ^-expansion is also expected to be 
powerful to analyze dependence on dimensionless parameters of solutions, like 
dependence on the number N and a ratio of two gauge couplings of U(N) = 
[17(1) x SU{N)]/Z n for an U(N) vortex. 

We expect that the ^-expansion can be applied to systems of non-critical 
coupling, although it might not be a straightforward extension. Our final goal 
is to establish a systematic tool to study the dynamics of vortices quantitatively 
without taking the critical coupling limit. Since in the ^-expansion vortices are 
treated as singular particles (strings) in a three(four)-dimensional spacetime, 
it will become possible to treat vortices of arbitrary shapes and discuss their 
dynamics analytically and quantitatively if we can consider such an extended 
^-expansion. 
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A Inequalities 

A.l Uniqueness of the solution 

Let us show the uniqueness of the solution f(x) of the following d-dimensional 
partial differential equation defined by a strictly increasing function W(/) with 
respect to / and a source term J as 

- d?f(x) + W(f(xj) = J{x), (A.l) 

where we require that /( x) vanishes at the spatial infinity. Note that if there 
exists a region T,f with its boundary <9£y for a certain scalar function f(x) so 
that f(x) satisfies 


f(x) < 0 for x £ £, f{x) = 0 for x € 3£, (A.2) 


which gives n-df(x ) > 0 with a normal vector n of c>£ and then Stokes’ theorem 
tells us the following inequality 



dS-df(x) > 0. 


(A.3) 


If we assume that there exist different two solutions fa (a?), fa(x) for Eq. (I A. Ill , 
then there exists the region £ 5 / for a difference Sf = fa — fa{or fa — fa) and 
we can derive inconsistency as, 

0 <f d d xd*6f(x)= f d d x{W(fa(x))-W(fa(x))} <0. (A.4) 

J'Ssf 


Therefore, if there exist a solution of Ea. (lA.l[l . then it must be unique. 

Furthermore, let us consider a solution f(x) with W(0) = 0 and J{x) > 0, 


-dff(x)+mf(x))> 0 . 


(A.5) 


If there exist a region £y for this function / where W(/) < 0, then we find 
inconsistency again 


0 < [ d d xd?f(x ) < [ d d xW{f{x)) < 0. (A. 6 ) 

J'Lf J'Lf 

Such a solution f(x) must be, therefore, positive semidefinite everywhere. 
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A.2 Sequence of sets of upper and lower bounds 

Here let us modify the inequality Eg. (12.371) for v > 0. 

: oo > xp > 0, 0 > P = r^- > —2v, (A.7) 

or 

to obtain a stronger set of upper and lower bounds of them. 

By integrating Taubes equation and P = rip' , we find relations between P 
and xp using integrals as, with Y = (r/Ri n ) 2v and setting m = 1, 


xp = ^[-P] = lira <j — 2 ^log ( jj— ) + I —P(s ) [■ = — logF 


e-s-0 


pr 

P = V[xp] = -2u + J dss (l — e“^ (s) J . 

Let us assume that the following set of inequalities T n 

Zn- f™>xp> /“, g™ > P > g™, for all r <E R> 0 - 
with some given functions satisfying 


• > /“ 1 > C > fn > fn—1 > > /o m =0, 

0 = g™ > ■ ■ ■ > Sri-l > g™ > 9n > 9n -1 > * * ' > 9o = “ 2l/ - 


2u), 


(A. 8 ) 


(A.9) 


(A.10) 


Using these inequalities, we can construct an another set of inequalities as 

*[ffn ] > xp > P[/f] > P > V\f™}. (A.11) 

Therefore we obtain a set of stronger lower and upper bounds as X n +i by 

9n+i = [tf, P[/ n M ]] , g ™ + 1 = max [g™, V[C }}, 

fn+i = min [/„ M , T[^]] , / n m +1 = max [/”, + [<?“]]. (A.12) 

Consistency of these inequalities requires that g^ + i > g ™+1 and f^_ 1 > f™ +1 
which reduce to, non-trivial inequalities 

0 = 5o M >TO, *[^]>/o m = 0. (A.13) 


This couple of inequalities turns out to give lower and upper bounds for R m as 
follows. 

The initial set of inequalities Iq gives 


Ii : oo > xp > max[ 0 , — logF], min 


— 2u, 0 


> P > - 2 v, (A.14) 


and therefore we find the followings are required 

f ' 

'0 


0 > ma xV[f* 


2v > 


drr 1 — 


r 

Pin 


2 l/' 


vRL 


2(1 + 1 /)’ 


0 < min 


Mi 


0 < —— log Y 


-=2 ^ 


= i/log 


R- e 
4 v 


(A.15) 
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and that is, R m must satisfy 


2y/is + 1 > R m > 2 \l — ■) 


(A.16) 


otherwise a function ip can not satisfy the set of inequalities Iq and thus blows 
up at large r. With R[ n satisfying the above set of inequalities, the next set of 
inequalities I 2 can be consistently obtained as 

f — log Y for r < ‘lyfv 

1 2 : m«[0 ,-logy]<^<| loir>2 ^ ■ 

(l-Tfc) 


r 2 

j 

[y-H 

> p> 1 


-2i/- 


-2v- 


vRf„ 


for r < Ri r , 
for r > Ri T , 


• (A.17) 


2 ( 1 + 1 /) 

In principle, you can calculate 13 , 14 ,..., sequentially as you like. 


B Some Integrals 


Since the modified Bessel function of the second kind is a two dimensional 
Green’s function, we can find the following relations 

J d 2 xKo(m\x — xi\)Ko(m\x — X 2 D 

27r 


— d 2 + m 2 


2 tt 


— d 2 + m 2 


A'o(m|fi - x 2 \) 


’ X 2t~ d 471-2 ^ 

a = Trf ' 5 


(9 7f 

= -27T-—-AA(m|fi - x 2 |) = — |xi - af 2 |A - i(m|fi - x 2 \) (B.l) 

om z m 


By using the integral formulas 


AT, 


° {x)= r^ e ~ Ht+i) ’ 


dO 

2 ^ e 


.tc cos 0 


one can calculate the following definite integrals, 


(B.2) 



’ r d 2 r 

drrI 0 (r)K 0 (r) 3 = J — e* 1 A' 0 (|£|)" 


d 2 a; dhdhdts cXl -(t 1+t2 +t 3 )- ti+dl (-L+-L+-L) 
2-7T 8tit2^3 

00 dtidt2dt^ 


0 4*1^3 


0 -(tl+t 2 +t 3 )+(^ + ^ + ^-) 


4 J 0 (1 + ui)(l + u 2 )(ui + u 2 ) 16 


(B.3) 
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with 1 1 = su\,t 2 = sii 2 ,tz = s , 



/ 


d 2 xd 2 y 

4w 2 


h 


(|f|)A' 0 (|f|)A'o(|f 


y\) K o(\y\) 2 


i r 

4 7o 

117T 2 

432 ’ 


dhdhdtadU e -fa + fa + « 3 +t4 ) +(^ + ^ RI7 ^_ =r ) 
HO + (U + 0)(f3 + ti) 


(B.4) 



0 


r d rrK 0 { r?= f_****** 

io 7o 4fH 2 f 3 (£ + £ + i) 


, — (tl+^2+^3) 


du\du2 


4 7 0 (1 + Ml + U 2 )(wi + M 2 + U 1 U 2 ) 


1 

36 


v-w (-) +v> (1) 


87T 2 

3 


0.585977 


(B.5) 


where ip^(x) = d 2 logT(x)/dx 2 is the digamma function, and 




drrA'o(r ) 4 


^C(3) ~ 1.051800. 
8 


(B.6) 
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